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Abstract 

We investigate a non-equilibrium ionization state and two-temperature structure of the intracluster 
medium in the linked region of Abell 399/401, using a series of N-body + SPH simulations, and find that 
there exist significant shock layers at the edge of the linked region, and that the ionization state of iron 
departs from the ionization equilibrium state at the shock layers and around the center of the linked region. 
As for the two-temperature structure, an obvious difference of temperature between electrons and ions is 
found in the edge of the linked regions. Ka line emissions of Fe xxiv and Fe xxv are not severely affected by 
the deviation from the ionization equilibrium state around the center of the linked region, suggesting that 
the detection of relatively high metallicity in this area cannot be ascribed to the non-equilibrium ionization 
state of the intracluster medium. On the other hand, the Ka emissions are significantly deviated from 
the equilibrium values at the shock layers, and the intensity ratio of Ka lines between Fexxiv-xxv and 
Fexxvi is found to be significantly altered from that in the ionization equilibrium state. 
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1. Introduction 

According to a standard scenario of hierarchical struc- 
ture formation in the universe, galaxy clusters are formed 
through successive merging of galaxies, galaxy groups, and 
clusters. About a half of observed clusters have irregular 
X-ray morphology (e.g., Akahori, Masai 2005 references 
therein), and a part of such irregularities thought to be 
caused by such merging events. 

Abell 399/401 is well-known as merging clusters on 
an early stage of the merging. Fujita et al. (1996) and 
Sakclliou, Ponman (2004) found that the intracluster 
medium (ICM) at the linked region of the two clusters 
is compressed based on the ASCA and XMM-Ncwton ob- 
servations, respectively. Recently, Fujita et al. (2008) re- 
ported the detection of Fe K emission lines near the center 
of the linked region based on the Suzaku XIS observation, 
and the metallicity in that regions is estimated to be 0.2 
times the solar metallicity, which seems relatively high 
compared with theoretical predictions based on numerical 
simulations (e.g., Tornatore et al. 2007). 

The estimation of physical properties of ICM in X-ray 
observations is usually based on the assumptions that ICM 
is in the ionization equilibrium state and that electrons 
and ions share the same thermal temperature. Such as- 
sumptions can be justified around central regions of galaxy 
clusters by the fact that ICM density is high enough to 
quickly achieve ionization equilibrium and thermal equili- 
bration between electrons and ions. As for the ionization 
equilibrium, the timescale required to reach collisional ion- 
ization equilibrium for an ionizing plasma, is estimated 



as n c t > 10 12 cm~ 3 s (see e.g., Masai 1984), where n e 
is the number density of electrons. Therefore, for n c ~ 
10~ 4 cm~ 3 in the linked region (Sakelliou, Ponman 2004; 
Fujita et al. 2008), t ~ Gyr is comparable to or longer than 
the merger timescale, so that the ionization equilibrium is 
no longer a reasonable assumption. Actually, in the warm- 
hot intergalactic medium (WHIM), where the density is 
much lower, the deviation from the ionization equilibrium 
have been pointed out by Yoshikawa, Sasaki (2006). In ad- 
dition, the thermal equilibration in merging clusters has 
been studied by Takizawa (1999) based on N-body/SPH 
simulations, and the two-temperature structure of ICM, 

1. e. the difference in temperature between electrons and 
ions, is reported especially at low-density regions. The 
two-temperature structure of WHIM has been also sug- 
gested by Yoshida et al. (2005). 

Therefore, there exists sufficient reasons to suspect that 
these assumptions are not valid in the linked region of 
Abell 399/401. In this paper, we investigate an ionization 
state and temperature structure of ICM in the linked re- 
gion by relaxing the assumptions of ionization equilibrium 
and thermal cquipartition between electrons and ions, and 
verify to what extent such a non-equilibrium state affect 
the interpretations of the observational data. 

2. Model and Calculation 

We carry out N-body/SPH simulations of two merg- 
ing galaxy clusters, in which a non-equilibrium ioniza- 
tion state and two-temperature structure of ICM are both 
taken into account. Radiative cooling and electron heat 
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conduction are both ignored in these simulations. 

The time evolution of two-temperature structure is fol- 
lowed with the same way as Takizawa (1999). We as- 
sume that the electrons and ions always reach Maxwcllian 
distributions with temperatures, T e and Tj, respectively, 
and the two temperatures are equalized through Coulomb 
scattering on a timescale of 
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where rij is the number density of ions, and In A is 
the Coulomb logarithm. By introducing the dimension- 
less temperatures of electrons and ions, T c = T c /T and 
Tj = Ti/T, respectively, normalized by the mean tempera- 
ture of electrons and ions, T= (n T c + niTJ) / (n e + n\) , the 
evolution of two-temperature structure is described by 
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where u and Q s h are the specific thermal energy and the 
shock heating rate per unit mass, respectively. The first 
term of the r.h.s. denotes the thermal relaxation rate 
between ions and electrons. Note that recent studies on 
supernova remnants showed that T c /Tj ~ 1 around the 
shock for relatively small Mach numbers (Ghavamian et 
al. 2007); electrons may be heated by some plasma waves. 
However, for a cluster merger, the Mach number of the 
shock is 'very' small (< 5). In this case, waves may not be 
generated and electrons may not be heated. In this letter, 
we ignore the effect of electron heating by plasma waves. 

The time evolution of ionization fractions of ions is com- 
puted by solving 



dt 
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k=l i=j+i 
~ Ujfj + CXj+lfj+l, 
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where j is the index of a particular ionization stage con- 
sidered, Z the atomic number, fj the ionization fraction 
of an ion j, Sij the ionization rate of an ion j with the 
ejection of i electrons, and oij is the recombination rate of 
an ion j. Ionization processes include collisional, Auger, 
charge-transfer, and photo- ionizations, and recombination 
processes are composed of radiative and dielectronic re- 
combinations. Ionization and recombination rates are cal- 
culated by utilizing the SPEX ver 1.10 software package 1 . 
Actual calculations are carried out in essentially the same 
way as Yoshikawa, Sasaki (2006), except that the reaction 
rates are computed using the electron temperature, T c , 
rather than the mean temperature, T, in order to incor- 
porate the effect of two-temperature structure. We solve 
the time evolution of each ionization fraction of H, He, C, 
N, O, Ne, Mg, Si, S, and Fe, but we focus only on iron. 

To reproduce the situation of Abell 399/401, the ini- 
tial condition of the two clusters is set up as follows. We 
consider a head-on merger (Oegerle, Hill 1994), and as- 
sume that their collision axis is perpendicular to the line 
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of sight. Each cluster has the same shape for simplicity 
since the differences of the two clusters in shape and size 
are slight (Sakelliou, Ponman 2004). We adopt the King 
profile for the initial dark matter distribution, in which the 
radial velocity dispersion, cr* r , satisfies a^ r = k^T v ^/ jj,m, 
where fee is the Boltzmann constant and T V j r is the virial 
temperature (the virial radius is r v j r = 3.0 Mpc and the 
cluster mass is M(r vir ) ~ 1.7 x 10 15 M ). For the initial 
ICM distribution, we first assume the temperature pro- 
file, T(r) = (T vir //?)cxp(-r/r vir ), with = 0.6, and the 
density profile is obtained by assuming hydrostatic equi- 
librium. The above parameters are chosen so that X-ray 
surface brightness satisfies the /3-model best-fit (Sakelliou, 
Ponman 2004), and that the spectroscopic-like tempera- 
ture, T s i (Mazzotta et al. 2004), reproduces the observed 
temperature at the unperturbed regions by Markevitch et 
al. (1998). The number of particles of each galaxy clus- 
ter is set to a half million each for dark matter and ICM. 
The initial relative velocity of the two clusters is set to 
1050 km/s to reproduce the observed temperature of the 
linked region in terms of T s \, and the initial separation 
is ~ 7 Mpc. As for the metallicity, a spatially uniform 
metallicity of 0.2 times the solar abundance is assumed. 
It should be noticed that the metallicity, Z, is very in- 
sensitive to the resulting ionization state of ions as long 
as Z <C 1 • We also assume T c — T\ = 1 and an ionization 
equilibrium state at the start of the simulation. 

3. Result 

In the rest of this paper, we present the results of a 
snapshot of the simulation in which the separation be- 
tween the centers of the two clusters is ~ 3 Mpc so 
that their configuration is consistent with that of Abell 
399/401. Figure la-le show physical properties on a 
cross section along the collision axis of the two clusters. 
The white-solid square with 1.5 Mpc on a side shows 
the region roughly corresponding to the observed area 
by Fujita et al. (2008). In this region, the local ICM 
density, n (= n c + ri;) ~ 3 x 10~ 4 cm -3 (figure la), and 
T$\ = 6.48 keV are in good agreement with the observa- 
tions (Sakelliou, Ponman 2004; Fujita et al. 2008). 

We find that the shock heating rate is typically a few 
percent of the adiabatic heating rate inside the linked re- 
gion. This means that both the electrons and ions are 
mainly heated by the adiabatic compression in almost the 
same rate with no significant shocks there. Accordingly, 
as can be seen in figure lc, the electron temperature is 
only a few percent lower than the mean temperature in 
this region. On the other hand, at the edge of the linked 
region (the white-dashed rectangles), there exist signifi- 
cant shock layers with mach number of 1.5-2. In these 
layers, we can see that the electron temperature is typi- 
cally 10-20 % lower than the mean temperature, and even 
lower toward outer regions. 

Figure Id and le show the ionization fraction of FeXXV 
and the ratio of the fraction relative to that in the ioniza- 
tion equilibrium state, respectively. Here, the ionization 
equilibrium state is calculated assuming that T c = T\. In 
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Fig. 1. Maps of (a) the ICM density in units of log cm -3 , (b) the mean temperature in keV, (c) the ratio of the electron temperature 
relative to the mean temperature, (d) the ionization fraction of Fe XXV, (e) the ratio of the Fe XXV fraction relative to that in the 
ionization equilibrium state, and (f) the ratio of the line intensity in 6.6—6.7 keV band integrated along the line of sight relative to 
that in the equilibrium state. The white-solid squares and white-dashed rectangles indicate the regions corresponding to the Suzaku 
observation by Fujita ct al. (2008) and the shock layers, respectively. 



the center of the linked region, the FeXXV fraction is typ- 
ically 30-60 %, the largest among other ionization states, 
and is 10-20 % larger than that in the equilibrium state. 
In the shock layers at the edge of the linked region, the 
FeXXV fraction is nearly 80 % and 30-40 % larger than 
that in the ionization equilibrium state. 

The excess of the Fexxv fraction can be understood 
as follows. In the linked region, the electron tempera- 
ture is increasing from ~ 3 keV to ~ 8 keV according 
to the adiabatic compression and/or the shock heating, 
so that the FeXXV fraction is decreasing in time because 
an Fexxv fraction peaks at ~ 3 keV and decreases as 
the temperature increases in the ionization equilibrium 
state. However, the ionization of Fexxv to higher ion- 
ization states is not quick enough to catch up with the 
ionization equilibrium state, leaving the FeXXV fraction 
larger than that in the equilibrium state. 

Deviations from the ionization equilibrium state can be 
seen also for Fexxiv. Its fraction is ~ 2-3 %, and ~ 15- 
20 % larger than that in the equilibrium state. Due to 
the excess of Fexxiv and Fexxv fractions, the fractions 
of Fexxvi and Fexxvn are ~ 5 % and ~ 15-20 % smaller 
than those in the equilibrium state, respectively. 

As a result of the deviation from the ionization equilib- 
rium, the intensities of iron line emissions are altered to 
some extent. We calculate X-ray spectrum from the sim- 
ulation using the SPEX package, and find that Ka lines 
of Fe XXIV and Fe XXV in a rest-frame energy band of 6.6- 
6.7 keV is larger than those in the ionization equilibrium 
state, primarily because of the excess of Fe XXV fraction. 
On the other hand, since the Fexxvi fraction is smaller 
than that in the equilibrium state, Ka lines of Fexxvi 



in 6.9-7.0 keV band in rest frame are dimmer. It should 
be noted that the iron line emission detected in Fujita 
ct al. (2008) corresponds to the Ka lines of Fe xxiv and 
Fexxv in the 6.6-6.7 keV energy band. Figure If depicts 
the ratio of the intensity in 6.6-6.7 keV band projected 
along the line of sight relative to that in the ionization 
equilibrium state, and it can be seen that the intensity 
changes are significant (typically ~ 15 %) at the shock 
layers, while the intensity is only a few percent enhanced 
around the area of the Suzaku observation by Fujita et 
al. (2008), despite the excess of the FeXXV fraction by 
10-20 % (figure le). This is because the deviation from 
the ionization equilibrium is significant only at the center 
of the linked region, and its effect is diluted in integrating 
along the line of sight. This suggests that the iron line 
emission detected by Fujita et al. (2008) is not severely 
affected by the deviation from the ionization equilibrium. 
Here, let us define the ratio of the X-ray intensity between 
the two energy bands as 



R = 



F(6.6-6.7 keV) 
F(6.9 - 7.0 keV) ' 



(4) 



Figure 2 shows the map of the ratio R/R cq , where R eq is 
the intensity ratio defined above but in the ionization equi- 
librium state, and clearly indicates that the intensity ratio 
departs from the one in the ionization equilibrium state 
within the shock layers. Note that this intensity ratio is 
independent of the local abundance of iron, but primarily 
depends on its ionization state. Therefore, such excess of 
the intensity ratio, R, can be used as a stringent observa- 
tional tracer of the non-equilibrium ionization state. 
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4. Conclusion and Discussion 

The ionization state of iron and two-temperature struc- 
ture in the linked region of Abell 399/401 is investigated 
by using N-body + SPH simulations by relaxing the as- 
sumptions of ionization equilibrium and thermal equili- 
bration between electrons and ions. 

It is found that, around the center of the linked region, 
the Fexxv fraction is 10-20 % larger than that in the ion- 
ization equilibrium state, and the electron temperature is 
only a few percent smaller than the mean temperature of 
electrons and ions, because the ICM is mainly adiabati- 
cally heated inside the linked region. On the other hand, 
we find that there exist shock layers at the edge of the 
linked region, and that the Fexxv fraction is larger than 
that in the ionization equilibrium state by 30-40 %, and 
the electron temperature is typically 10-20 % lower than 
the mean temperature around these layers. 

Our simulation indicates that the intensity of Fe Ka 
emission lines is affected by such deviation from the ioniza- 
tion equilibrium state of iron in the linked region. While 
the deviation from the ionization equilibrium state is re- 
markable around the center and the edge of the linked 
region, we find that the emission line intensity is strongly 
affected preferentially at the edge of the linked region. 

The area of the Suzaku XIS observation by Fujita et 
al. (2008), in which fairly high metallicity is reported, is 
located at the central portion of the linked region. Our re- 
sults imply that X-ray emission in this area is not strongly 
affected by the effects of a non-equilibrium ionization state 
and two-temperature structure. Therefore, the fact that 
high metallicity is detected in this area cannot be ascribed 
to the non-equilibrium ionization state of the ICM, and 
must be explained by other physical processes. 

It is interesting to discuss the detectability of the 
non-equilibrium ionization state and the two-temperature 
structure of ICM in Abell 399/401. According to our sim- 
ulations, it is expected that shock layers with a mach 
number of 1.5-2 are located at the edge of the linked 
region, and that the ratio between Ka emission lines of 



Fexxiv-xxv and FeXXVl are significantly different from 
that in the ionization equilibrium state. Observationally, 
such an intensity ratio of iron emission lines could be 
important clues for the detection of the deviation from 
the ionization equilibrium in these layers. In principle, 
the indication of the two-temperature structure can be 
also obtained by the difference between electron and ion 
temperatures inferred from X-ray thermal continuum and 
thermal width of emission lines, respectively. Of course, 
the detections of the non-equilibrium ionization state and 
the two-temperature structure are not very feasible with 
the current observational facilities, but could be achieved 
by the X-ray spectroscopy with a high energy resolution 
using X-ray calorimeters in near future. 

Finally, we should discuss several caveats on the as- 
sumptions adopted in this work. The slope of the den- 
sity profile, i.e. /?, is a sensitive parameter to the non- 
equilibrium ionization state and two-temperature struc- 
ture. If we adopt [3 = 0.5, the resultant departure from 
the ionization equilibrium is suppressed because the ICM 
density at the outskirts is relatively denser and the equili- 
bration timescales become shorter. In addition, we assume 
that the encounter is taking place on the plane of the sky. 
If it is not the case, the effect of the non-equilibrium ion- 
ization state on the observed X-ray intensity at the shock 
layers would be blurred to some extent in integrating along 
the line of sight, depending on the viewing angle. 
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